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'nI" ■ I. INTRODUCTION 

<^ . 

Nickel [l| found singularities in the susceptibility of the two dimensional isotropic Ising lattice in the absense 
of an external magnetic field. Explicitly, he found that for the isotropic Ising lattice with lattice constant K, the 
susceptibility has a natural boundary on the circle 



sinh2ii'| = l, (1) 



unlike the free energy and the magnetization. That is, the set of singularities of the susceptibility is dense on the circle 
'~~'> ([1]), whereas the sets of singularities of the free energy and the magnetization are not; in fact, they are dense only 
' at isolated points. The reason for this is not yet known, but there are som suggestive findings regarding differences 
^ between the susceptibility on the one hand and the free energy and the magnetization on the other. A dijferentiably 
[■ — finite function is one that is the solution of a linear ordinary differential equation of finite order with polynomial 
0^ , coefficients. Guttmann and Enting Q found that while (in the absence of an external magnetic field) the partition 
function and the magnetization are both differentiably finite functions, there is strong evidence that the susceptibility 
is not. They also found that the susceptibility most likely has a natural boundary even in the non-isotropic case. 
■ Specifically, what they found was that for an anisotropic lattice with lattice constants Ki and K2, the susceptibility 

O ' x(-2i: Z2), where 
00 . 

O' := tanhiTi (/ = 1,2), (2) 

• ^ , has a natural on the circle \zi\ = 1 for fixed Z2- Further evidence supporting this claim was given by Orrick, Nickel, 
^ ' Guttmann and Perk [1]. 

^ ] The Haniiltonian of the Ising lattice of size (M, N) in an external magnetic field H is 

_ _ I 

M N M+1 N M N 

£ ■■= -El ^ ^ 0-j, fecjj- fc+1 - ^2 ^ ^ (Jj,k(Tj+iM -H^^ aj.k, (3) 

j=l k=l j=l k=l j=l k=l 

where no boundary conditions have been specified. Denote by Zm,n the corresponding partition function. Exact 
solutions may be found for the particular magnetic fields H = and H/kT — i-n /2. In what follows, the dimensionless 
parameters Ki :— Ei/kT {I — 1, 2) and h :— H/kT will be used instead of Ei, E2 and H. It is customary to use the 
variables 

a;; :=e-2^' (Z = l,2) (4) 

and 

z (5) 

instead of Ki, K2 and h. 



* Electronic address: 'i var.lyberg@uclouvain.be] 



2 



It was conjectured by Fisher [s] that in the isotropic case K Ki ~ K2, the zeros of Zm,n in a; := e wiU all 
approach the circles 

\x±l\^V2, (6) 

or equivalently the circle ([T]), in the thermodynamic limit AI, N ^ 00. Lee and Yang showed that if < < 1, then 
zeros of Zm,n in z all lie on the unit circle \z\ = 1. Wood 7] showed that in the non-isotropic case {Ki, K2) = {K, aK) 
(a S R) , the zeros of the partition function will all approach the curve described by the polar equation 

- 11 + r2 

cos 6* = — ^ — (7) 

in the thermodynamic limit, where 6 is the polar angle. In particular, this means that the zeroes of the partition 
function in the anisotropic case with lattice constants Ki and K2 will spread over a two dimensional region (as a is 
varied). It should thus be possible to analytically continue the free energy through this region. The transformation 
K2 — > —K2 reflects the curve ([7]) about the imaginary x axis. Thus, in particular, for a = 1 these two curves are the 
two circles 

In the thermodynamic limit Ad, N 00, the free energy 

f{z,xi,X2):^kT lim --^\ogZM,N (8) 

will only depend on z, xi and X2. Away from the zeros of the partition function, / will be an analytic function of xi, 
X2 and z. 



A. The Brascamp-Kunz lattice 

The partition function Zm.n of a finite lattice with the Hamiltonian ([3]) has been calculated exactly for many 
different boundary conditions, including for the Ising lattice on a Mobius strip and on a Klein bottle Q. However, 
the boundary conditions particularly useful for studying the natural boundary are the so called Brascamp and Kunz 
boundary conditions, described in section [ill Brascamp and Kunz [9], using a result by McCoy and Wu ^M^, showed 
that for the Ising lattice with so called Brascamp-Kunz boundary conditions, the zeros do not just approach the set 
([6]) asymptotically, but lie precisely in it. Brascamp and Kunz constructed the boundary conditions named after them 
from McCoy's and Wu's calculation of the partition function of a cylindrical lattice with a magnetic field on the 
lower boundary. McCoy and Wu considered a general ^ (although strictly speaking Sj/kT ^ i7r/2), but Brascamp 
and Kunz found that in the limit f)/kT — > J7r/2, the expression for the partition function much simplifies, and the 
zeros of the partition function all lie on the circle ([T]). They then used the duality relations presented in ref. [ll[ to 
construct a lattice whose dual lattice is the lattice of McCoy and Wu with f)/kT = 111/2. This lattice is the Brascamp 
and Kunz lattice. Specifically, it is the two dimensional cylindrical Ising lattice A = {(m^n) | 1 < m < M, 1 < n < 
N, (to, + 1) = (to, 1)} with the following boundary conditions. 

(i) The lattice interacts with a row of fixed, positive spins above it, 

(ii) The lattice interacts with a row of fixed, alternating spins below it. 

(See fig. d]). The interaction between neighboring spins on A is Ei in the horizontal direction and E2 in the vertical 
direction. Apart from from interactions between nearest neighbors, there is an external magnetic field H{j, k). In this 
paper, H{j, k) will constant; either H{j, k) =Q for all (j, k) or H{j, k)/kT — «7r/2 for all (j, k). Thus the Hamiltonian 
is defined as 

M N M N M N 

fA(cr, Ei,E2, H) = ~Ei ^ ^ crj,fcCrj,fe+i ^ -^2 ^ ^ f7j,fcCrj+i,fe - X! X! ^(■?' (9) 

j = l k=l j=0 k=l j=l k=l 

where Uj^k = (^],k+N = ±1, cro,fc = 1 and aM+i.k (-1)*'+^ 



B. The magnetic field H/kT = i-K/2 



The two dimensional isotropic Ising lattice in the presence of an external magnetic field H/kT = z7r/2 was first 
studied by Lee and Yang Since then, this problem has been much studied, both for isotropic [13 [Si and 
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FIG. 1: The lattice A (marked by "x"; M = 3, A'' = 4) and the dual lattice A* (at intersections of lines). 

anisotropic lattices[l3|. However, to the authour's knowledge, no study has been made of a lattice with Brascamp- 
Kunz boundary conditions in an external magnetic field. 

As will be shown in section |TT1 the set of zeros of an isotropic lattice with these boundary conditions in a magnetic 
field H/kT — iir/i will lie on the one dimensional, connected set; moreover, in the thermodynamic limit the set of 
zeros will become dense on this set. In this sense the magnetic fields H = and H/kT = i7r/2 are similar. It has 
been conjectured elsewhere [31 that for any other magnetic field, the set of zeros does not become dense on a set that 
has even a one dimensional subset. Future studies may reveal the relationship between Nickel singularities and the 
particular magnetic field H = 0. It is yet unclear whether there are Nickel singularities for iJ ^ 0, or only for H ~ 0. 

Lee and Yang [S] showed that the isotropic Ising lattice in a magnetic field z7r/2 is equivalent to a two dimensional 
lattice gas. Let the free energy per site in the thermodynamic limit be /. Lee and Yang conjectured that the 
magnetization I{i7r/2, K, K) := {kT)~^{df /dh)\i^^i^/2 is 



I{iTT/2,K, K) 



where 



1 — a;^ 



-2K 



1/4 



(10) 



(11) 



Later McCoy and Wu [T^ proved (fTOj) by calculating hm„^oo •S'n =: P where the row correlation function 

Sn '■= (co.oco.n) IS the expectation value of the product cto.oCo,™ of spins at sites (0,0) and (0, n). Moreover, they 
calculated the magnetization for the non-isotropic lattice with lattice constants Ki and K2 in the horizontal and 
vertical directions, respectively. They found 



I{i7r/2,Ki,K2 



1/4 



(12) 



where 



zi 



:= tanhif;, 1 = 1,2. 



(13) 



Gaaff [I4I found that the Ising lattice with an external magnetic field H is in fact equivalent to a sixteen vertex model, 
which simplifies to an eig ht vertex model when H/kT = m/2. Wu [11], using the same approach as Gaaff, found 
that the Ising lattice has no phase transition in a magnetic field H/kT = i7r/2, unless additional diagonal interactions 
and/or four-spin interactions are introduced. 

While numerical studies have been made on the problem of the Ising lattice in a general magnetic field [3 ■, there are 
no theorems on the analyticity of the free energy in this case. It should be possible to obtain the analytic properties 
of the free energy / from the zeroes of the partition function Zm,n- f can only depend on Ki, K2 and h. When 
Ki = K2 = K , one thus seeks the zeroes of Zm,n iii the planes of z and x respectively. However, no systematic study 
of the zeros of the partition function has yet been made. The finding that the set of zeros becomes dense on a one 
dimensional set both ioi H = and H/kT = in/ 2 may be useful for such a study. 
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II. SUMMARY OF RESULTS 



The partition function 



Z^{K^,K2,h)= e,^p^£A{a,Ei,E2,H)/kT (14) 

CTe{-i,i}A 



and 



will be calculated for the constant external magnetic fields H = and H/kT = z7r/2 for which the problem is exactly 
solvable. Brascamp and Kunz calculated Za{K,K,Q). Za{Ki,K2,0) and Z\{Ki, K2,in/2) are given by 

N/2 M 

Za{Ki,K2,0) = x;^'^^\-^'^^^ n + + ^2) - 2x2(1 - xDcosOj - 2x1(1 - xDcos^k} (15) 

j=l k=l 

ZA{K,,K2,^./2) = u-''^'\2'''"' 

N/2M/2 

n ri'f^^ + + "2) -4uiU2 - 2u2(l -ui)^cos20j - 2ui(l -U2)^cos0fc}, (16) 

3 = 1 k=l 

where 

Oj = (2j - l)^/N, ifk = k^/{M + 1), 0fe = (2fc - 1)^/(M + 1) 
and ui = e-'^^' {I = 1, 2). (17) 

In particular the free energy per site in the thermodynamic limit, /, is 

log{(l + ul)(l + ul) - 4uiU2 + 2wi(l - U2)^ + 2^2(1 - ui)^ - 4^2(1 - wi)^cos^x - 4ui(l - U2)^ cos^ y}, (18) 

which is the same as the result found by McCoy and Wu 15]. 

As has been already mentioned, with Brascamp-Kunz boundary conditions, the zeroes of Z/^{K,K,0) lie on the 
circle even when the lattice is finite, as can be seen from p5|) . In the case h = «7r/2, the set of zeroes will also 
be one dimensional with Brascamp-Kunz boundary conditions in the special case Ki — K2. When Ki — K2, (fTB]) 
becomes 

N/2 M/2 

Za{K, K, ^7^/2) = {u-^ - 1)^^/^ H IT + + " ^^^^^ 't"'/^ ~ ^ ^j)^ (1^) 

3 = 1 k=l 

The zeroes lie on the circle 

|u| = 1 for - 1 < 3 - 2 cos^ Oj - 2 cos^ (t>k < 1 (20) 

and on the line segment 

- 3 - 2V2 < u < -3 + 2V2 for 1 < 3 - 2 cos^ 0j - 2 cos^ (pk < 3. (21) 



III. THE CALCULATION OF Za{Ki,K2,0) 
It will be shown that the partition function Za{Ki, K2, 0) is given by p^ . 

Brascamp and Kunz '01 used a result found by McCoy and Wu [13] to get their result. Ref. [13] concerns a 
lattice on a finite cylinder with a magnetic field on the lower boundary. To prove (|15|). we consider the dual lattice 
A* = {(m, n)\l<m<M + l, 1 < n < N, (m, N + 1) = (m, 1)}, shown in fig. [T] It is periodic in the horizontal 
direction and free along the upper boundary. The lattice constants are K2 in the horizontal direction and Ki in the 
vertical direction, where 



sinh 2Ki sinh 2i^* = sinh 2K2 sinh 2K; = 1. 



(22) 
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Therefore there may be defined a Hamiltonian on A* given by 

M N M+l N M+1 N 

(a, Kl,K*,H*) := -E*, ^ ^ tr,- - i?2 E E '^o^^^J-M+i - ^*(^' (23) 

j=i fe=i j=i fc=i i=i fc=i 

and a corresponding partition function of A* given by 

ZA.{KlK;,h*)^ E exp-£A.(a,A7,X*,i7*)/fcr*. (24) 
cre{-ia}A* 

In order to continue, it is necessary to establish the following lemma: 
Lemma 

ZA(ifi,X2,0) = 2-^-''/\smh2K^Y''''\sixv\i2K2)^''+^^^'^ZA.{Kl,K;, [)) (25) 
where \) is an external magnetic field on A* given by 

\){ni,n) = —5m.M+i- (26) 

Consider fig. [2] There is a perpendicular line drawn between any two spin sites on A whose spins have product 
— 1. With this particular spin configuration, each site on the lower boundary of the dual lattice A* has precisely one 
dimer intersecting it. It is easy to see that with any spin configuration, each site on the lower boundary of A* will 
have one or three dimers intersecting it. Let a dimer carry the weight of the lattice constant on A that it crosses, and 
let the collection of dimers corresponding to a spin configuration be called A. The energy of a general configuration is 

£;(A) = -MNEi - {M + l)NE2 + 2|A| (27) 

where |A| :— hE2 + vEi and h and v are the number of horizontal and vertical dimers, respectively. In this way the 
partition function may be written as a "low temperature" expansion 

ZA(ifi,X2,0) = ^e-2|>IAr, (28) 

A 

where the sum is over all possible paths A. 

Now consider the dual lattice A* on its own. A* is a cylindrical lattice with free boundary conditions. Suppose the 
magnetic field [} is applied. Since exp(Ti7r/2 = icr, it follows from ([^5]) and ([M]) that 



Za. {ki,k;, [)) = E ( n J cxp -8a. (a, kik;,o). 

Since exp K^aa' — coshifj* + crcr' sinh isTj* , this can be written as a "high temperature" expansion 
Za.{KI,K;,1)) = (cosh if*)*^^ (cosh X*)^^^+'^^ 

CN \ / M N \ /M+l N \ 

fe=l / \j = lk=l I \m=l n=l / 



(29) 



(30) 



TG{-1,1} 

where z* := tanhif*. Rewrite the sum as 



where 



Za-^KI^K^^'c)) = (coshif*)*'^^(coshX*)^^+^^^ 
/w/2 \ / 

E n «^220'M+l,2fe-lO-M+l,2fe (1 + K(a;,2;)crj;Crj,) ] , (31) 

<Te{-l,l}A* / \{a;,aeA- I \x-y\ = l} 



zl if - = (±1,0), 

n{xj,yj) := \ z* if - = (0, ±1) and ixj,yj) ^ ((M + 1, 2fc - 1), (M + 1, 2k)) for each fc, (32) 
l/z^ if Xj - yj = (0, ±1) and ixj,y/) = {{M + 1, 2fc - 1), (M + 1, 2k)) for some k. 
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FIG. 2: A spin configuration of tlie lattice A witfi no external magnetic field. Here M — 3 and A'^ = 4. A perpendicular line is 
drawn between any two spins of opposite signs [TH . In this way, this particular configuration gives rise to one closed loop and 
two dimers on the dual lattice A*. The latter give the appearance of a magnetic field S)* /kT* = i7r/2 on the lower boundary 
of A*. 

Clearly, the only terms in the sum that contribute to Za» {Ki, f)) are of the form 

Z2'^M+l,iaM+l,2Z2(^M+l,3CrM+l,4---Z2CrM+l,N-lCrM+l,N 

[K{x„+i,yn+l)(J.j:^ + iay^^^K{Xn+2,yn+2)(^x„ + 2^y„+2--l^{Xp,yp)(Txp(Typ] 
■■\l^{Xq+l,yq+l)(J.^^^^ay^^^K{Xq+2,yq+2)(yx,+2<^y,^2---l^{Xr,yr)<7x^Cry^] 

or else of the form 



^2 0'M+l,20'Af+l,3^2 0'Af+l,40'A/+l,5---^2 ^M+l.N^^M+l,! 
[K(xi,yi)a.j;,(Jy^K{x2,y2)crx2<^y2--'^iXn,yn)<J 

[KiXn+l,yn+l)(7x„ + i(Ty^ + iK{Xn+2,yn+2)<Jx„ + 2^y„+2--K{Xp,yp)ax^ayJ 
...[K{Xq+l,yq+l)a.^.^^^ay^^^K{Xq+2,yq+2)(Tx,+2(^y,+2--l^iXr,yr)crx^(Ty^] 

where Xj and yj are nearest neighbors, yj — xj+i and xi — yn, Xn+i = yp,---, Xq+i = yr- Let each pair {xj,yj) be 
marked by a dimer of weight K{xj,yj). Should two dimers have the same position, then clearly one will have weight 
^2 and one l/z2, and thus they will form one dimer of weight 1; that is, no dimer at all. Let a contributing path 
as above be called 7, and let ||7|| := {zl)'"{z2)'^ where v is the number of vertical dimers and h is the number of 
horizontal dimers. It follows that Za*{K*,K2, ()) can be written as 



Za. {K*,K;, t)) - 2(*^+i)^(coshif *)*''^(coshi^*)^*'+'^'^ Il^l 

7 



where the sum is over all contributing paths 7. Since tanhX;* 



, it follows that 



^"7|| = 2Ve-'l^l/'=^ 
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A 



(33) 



(34) 



and hence the lemma follows. 

Remark The general argument concerning duality relations was first given in ref. 

McCoy and Wu did not restrict the magnetic field on the boundary, i^, to any particular value, so the calculation 
done in ref. \ld^ is more general than necessary to solve our particular problem. What follows is a proof of (jl5[) . 
It follows from © and HH) that Za{Ki, K2,0) is 



ZA(ifi,X2,0) = (cosh Xi)*^^ (cosh ^2)^*'+'^ 



M N 



(yj,k<Jj+l,kZ2) ) ]J(1 + <Jl,lZ2) 



N 



(35) 



1=1 
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where z/ = tanh Ki . In the same way it can be shown TT| that 
ZA'iKl,K;,i)) = (coshA7)*^^(coshX*)*^''+^^^ 

X! ( II n^-^ + ) ( II n^^ +'^'"■"'^'"^"+1^2)) 11*^*^+1^'- ^^^^ 

o-e{-l,l}A* ^j = lfc=l ^ \ m=l ri=l ^ i=l 

It remains to show that ([TS]) holds. The calculation of the partition function Za» {K*, K2, f)) is just a special case 
of a calculation already done by McCoy and Wu 10]. Since sinhi7r/2 = i, it follows from equation (3.26) of ref. [lol | 
that 

Za. iK;,K;, [)) = 2(^^+i)^(coshif*)(*^+i)^(coshifr)*^^Zi**^^/'(l - )(a^+i)a^/2 
II 



2(^^+i)^(cosh/f*)(*^+i)^(coshif*)*^^z*^'^^/2(l - zf 

n "jP^K--«r') - 2cos(pfc +a,-i), (37) 



where aj is the larger root of the quadratic equation 

(1 + zf){l + zf) - 2z^(l - zf) COS0J - zl{l - zf){a + a-^) = 0. (38) 

Explicitly 

aj = 2- {(1 + zf){l + zf) - 2z*(l - zf )cose, 

+ (1 - zf )[{\ - aie^^O(l - aie-''>^){\ - &^^e'<'')(l - a^^e-^'^O]'^'}, (39) 

where 

«1 = Z*(l-Zi*)/(1 + Zi*), 

«2 = zl-\l~zl)l{l + zl). (40) 

Oj and v'j are given by 

_ n A |l + z,V'^^p-4zf sin^g,|l + z,V^^|-2-(l-zf)2|l + z,V'^^|-2 \y/^ 

^ V2 V (l-zf)[(l-die'»0(l-aie-*«0(l-a2"^e»''0(l-a2"^e-'«0]^/^ 

and 



Since 



^, ^ / 1 A ^ zf |1 + zle^<^^ p - 4zf sin^ e,\l + z*,e^<^^\-' - (1 - zf )2|1 + z^e^<^^\-' \ \ ^ ^^^^ 



(1 - zf )[(1 - aie'»0(l - aie-^«0(l + a2"^e*»0(l ' aa^^e-^^Ol^^^ 



1-zf 



z*(l-z* ) 
it follows from ([57]) that 



-1 ^*.X JV/2 

Za.(Xi*,X*,[)) = 2(*'^+i)^(coshif*)'*'+'''^(coshXi*)*^^ ' 1 
2sin^ 



^ bill (/t 

M 

H {(l + zf)(l + Z2*^)-2z2*(l-zf)cos0j-2z*(l-Z2*^)cos(^fc}. (44) 



fc=i 



Since 

cosh 2if;7 cosh 2is:, = smh2A';, 
|l + e' 



M , ^^e,^2 
2sin^ 



J = l 
N/2 

and n - " Aie-^^0(1 " e^'Va2)(l - e-^^Vda)]^/' - 2(1 ~ zfy'\zlz*)'''^, 

it follows from (jU]) that 

ZA.(i^*,if*,f)) = 2(*^+i)^+i(coshif*)('^^+'^'^(coshi^*)'^'^4'^^' 

N/2 M 

n n^(l + + ^2^) - 2^2(1 - 4^) COS 6lj- - 2z*(l - zf)cosipk} 

3 = 1 k=l 

= 2*^^+^/2+i(sinh2i^i*)^-^^/2(gi^j^2i^*)(A/+i)iv/2 

Jl {cosh 2i^i cosh 2K2 - sinh 2/^ cos 9j — sinh 2iir2 cos t^fc}. 

i=i fc=i 

It follows from (121, (1251) and (Si) that (HI) holds. 



IV. THE CALCULATION OF Zf,{Ki, K2,iTT /2) 
The partition function in this case can be found from ^ and lfT4)) to be 

ZA{Ki,K2,i7r/2) = (coshifi)*^^(coshif2)^*^+'^^ 

. M N N 
( n n + 0-j,fcCrj,fc+lZi)(l + 0-_,-,feCTj + l,fcZ2) j ]^(1 + (Tl,iZ2) 

where zi = tanh ii'; . This can be written as 

ZA{Ki,K2,i7T/2) = (coshifi)^^^(coshif2)(*'+')'^4''^/' 

, M N N 



where 

£0") ^ I ^2 if ^'^'5^' 
2 [ l/z2 if j is odd. 

(|49)) gives interactions Ki, K2 and on A, with the definition of K2 being 

tanh K2 = 1/ tanh K2 ■ 

This defines a new Hamiltonian 

M N M/2 N 

^Aio-, El, E2, E2, H) = --El ^ ^ 0-j,feO'j,fc+l - -^2 ^ ^ 0'2/,fcCr2/+l,fc 

j = l k=l 1=0 k=l 

A//2-1 N M N 

X! X!^2/+l,fcO-2/+2,fc fc)crj,fc, 
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FIG. 3: The interactions on the dual lattice A*. Here M — N — 4. This figure shows the three kinds of bond. The bond 
Ki is shown as a broken line, K2 is shown as a solid line, and — JfJ is shown as a dotted line. There is a magnetic field 
Sj* /kT* = in/2 on the lower boundary. 



and a new partition function 



<je{-i,i}A- 



In this way, can be written as 



fsinhiiT 
(cosniir2) ^ 



The Hamiltonian on A* is 



M+l N M/2 N 

j=l k=l 1=0 fc=l 

Af/2-1 N M+l N 

/=0 k=l j = l k=l 

and its corresponding partition function may be defined as 

Za- (if* , K; ,K;,h*)^ Y exp -£a' (a, El , £;2* , E; ,H*)/kT\ 
(Te{-i,i}A* 

where A* is the dual lattice (see fig. [3]) and 

sinh 2X1 sinh 2if 1* = sinh 2i!:2 sinh 2/^2* = sinh 2/^2 sinh 2^2* = 1- 
In the same way as (j25p was obtained, one obtains the relation 

Zj,{K,,K2,k2,Q) = 2^^^-i (cosh Xi)^^^ (cosh X2)(*^+i)^ 



[(tanhXi)*^^(tanhif2)('''/'+'''^(tanhi:2)*''^/']'/'iA* (i^i*, ^2*, ^2*, t))- 



ft follows from ^1 and IMI) that 



tanh K2 = — tanh K2 ■ 
Therefore the lattice constant K2 on A* may be taken to be such that 

K*2 = -if*- 



(52) 



(53) 



(54) 
(55) 
(56) 

(57) 
(58) 
(59) 
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Hence ([57)1 simplifies to become 

iA(i^i,if2,i^2,0) = 2^^^/2-^/2-l(sinh2Xl)*^^/2(sinh27^2)'^/'(coshi^2)^'^^^A.(i^l^i^2^i^^^ (60) 
ZA*{Kl,Ki,k^,t,) is given by 



AT 



ZA.{K*,K;,K;,i)) = (coshif*)*^^(coshif*)^*^+'^^ 

, M N , M+l N -.11 

( JJ(1 + crj,fcfTj+i,fcZ*) j f ]^(1 + 0-m,„crTO,„+iZ2*-'"'') j JJif7Af+l,i 

crg{-l4}A' ^j=lk=l ^ m=l n=l ^ 1=1 

where z*^"^ = (-l)™+iz2*. Let e > 0, and let 

= i7r/2 + e and = tanh/ig. 



The product n|!=i i(^M+i,i can be written as 



N 



N 



W_iaM+ij = lhn(cosh/ie)^]^(l + aM+i.iz^)- 



1=1 



1=1 



As shown in ref. y/Tj , the sum ^ of (pTjl can be written as 

^ ^ ^{M+l)N^l lini(cosh/i,)^PfA, 



(61) 



(62) 



(63) 



(64) 



where is a matrix which will be given below, and Pf^e is the Pfaffian of A^. The determinant is given by [13] [3l 



N 



del A, = J|detQ5'^)(6ij ) 



where S^'K^j) 

is a 4(Af + 2) X 4(Af + 2) matrix given by 



^Tn^m(^j) 



/ 

-1 - (-l)™+iz|e-'''^ 
1 
1 



1 + (-l)™+iz*e'''j -1 -1 \ 



1 -1 
1 
-1 j 



for 1 < m < M + 1, 



( 1 + e*^^ -1 -1 \ 

1 -1 

1 -10 1 

V 1 1 -1 y 



(65) 



(66) 



(67) 



for 1 < m < Af, and 



03 



^m+l,m(^j) 



/O 

. _ 

~ 

\ 

/O 





Vo 



All other elements are zero. Let 



a± — ±2iz2 sin 6*^11 ± Z2e 

b± = ±(l-Z2*')|l±Z2*e*' 

c = 2isin6lj|l + e*^^|-2. 



*„i6ii 1-2 



(68) 



(69) 



(70) 
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Using the same procedure as in ref . [10| , one finds that det A is given by 

JV 

det A, = [| 1 1 + e'^^ n 1 + z* e'^^ I 1 1 - e'^^ I det C^^^ (0j ) , 



(71) 



where 



/ -c 

C Ze 

-Ze -a+ b+ 

-b+ a+ zl 

-z\ -a_ b_ 

-b_ a_ z\ 



\ 



-b+ a+ / 



(72) 



The dimension of C^'^ is 2(A/ + 2) x 2(M + 2). Let <6n^ be the determinant of the 2(n + 1) x 2(n + 1) matrix 
form (|72p . and let J)^'^'' be the (2n + 1) x (2n + 1) determinant with the last row and the last column removed. 



^l-'-'2m+l 



and 



where 



b+ a+zi* 



-al + hi 



a-Zi 

_*2 



*-'2m-l 

-'I -^2)71-1 



1 < m < M/2 



zl^\ 



*2 
'-1 



Z^ 



Co 

,2,*-l 



Co = -c and 



74t becomes 



The matrix 



-a^ + b2 



„*2 



-al + b^ 
-a_z* 



0(1) 
0(1) 



of the 
Then 

(73) 



(74) 



(75) 



(76) 



a-zf 

*2 

'1 



^1 



1 - 4zi*''zr sin^ 6*^11 - zre'"^*^ 2izi*(l - Zi*^)z^ sin6lj|l - z^e''^^ \ 
2iz*(l — 21^)^2 sin^^jjl + Z2e 



iffj 1-2 



-2 ,*4 



4zJ'^Z2^ sin^ ( 



1=^ .1 
. — Zo e 



(77) 



has eigenvalues 



A = Zl (1 - Z2 ) a/|l - Z2 e ^ | 



and A' = zf (l-Z2*')Va|l-Z2*'e'2«^f 
where a is the larger root of the quadratic equation 

(1 + z*'')(l + Z2*^) - 4zf Z2*^ - 2z2*^(l - zf )2cos26»j - zfil - zff{a + l/a) = 0. 
*P has corresponding eigenvectors vi and W2- The eigenvector with eigenvalue A is 



(78) 
(79) 



and the eigenvector with eigenvalue A' is 

1 

It follows from JTSD, (HH) and ([771) that 



^2=1 62 ' " 



(^) \ 

1 = -z^ 



Since 



it follows that 



M/2 _ f ai a2 W A*^/2 \ / ai 02 



61 62 M M ^1 



Therefore 



Since 



CVi \ .2 / ai a2 \ / A^V2 Waia2\-7a+\ /o(l) 
^rsiti j " ^ I ^1 A A'^^/2 )[b^ h) [zl ) + [ o(l) 



2 

dete:(^) = e:^';^, , = ^ {aiA^/^(b2a+ - as^i*) + a2A'*^/2(_^^„^ 

ai02 - 0201 



a2(-&ia+ + fli^ij') a ^01(620+ - 022:1), 



can be written as 



dete:(^) = - ■ ai(62a+-a2.r) ;f' ^^2 "''^'("^''^''^' + "''''^'^^') + "(1) 

aib2 - 0201 |1 - Z2 e«2»j \M 

01(620+ - 02^1*) ,' ,2 Je- + n + - 2cos0fc) + 0(1). 



01&2-0261 |l-z^^e I 

Using (|86p again, one obtains 

ai(62a+ — 022:2)(1 + a^"^) = a+. 



01^2 - 02^1 

Hence ([57)) can be written as 

*M/-, _ *2nM ^'^/2 

det = -^^a+ ; ^ 4 ^,^ n + 2 cos<^.) + o(l). 

I-L ^2 e ^1 j,^^ 

([7T| and jlSI) give 

N M/2 

det A, = ]J|l + z*e'^^f |l + e*«:'|2z2ca+^*M^^^^*2^Af -Q^^^^-i _2cos0fe) + o(l). 
j=i k=i 

Using jTO]) and ^ can be written as 

N 

detA^ = Jj4z2sin^6'jZ^^ 
i=i 

M/2 

{(1 + zl^){l + zt) - 4zfzf - 2zf{l - zff COS20J - 2zf{l - zff cos^fc} + o(l). 
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(HZl), (EH), (El and give 

N/2 

Za,{K*,K*,K;,1)) = {cosh K*f^^^ {cosh K;)^^^+^^^Y[4siB:'9jz; 

M/2 

'[[{{l + uj){l + ul) - AuiU2 - 2?i2(l - uifcos29j - 2ui(l - itz)^ cos (92) 



fc=i 



Since 



N/2 

]j2sin6'j=2, (93) 



it follows that 



ZA'{KlK;,k*,i)) = 4(coshi^i*)*^^(cosliif*)'*^+'^^^2^^' 

N/2 M/2 

n {(1 + u?) (1 + u^) - 4uiU2 - 2u2(l - uiY COS 29 j - 2wi(l - cos 0fc}. (94) 
j=i fe=i 



It now follows from ^ and ^ that ([TBI) holds. 

V. CONCLUSION 

The new results ([T5|) and (fT6|) have been presented. 
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